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It is shown theoretically that a persistent current can be continuously created in a Bose-Einstein 
condensate (BEC) of alkali atoms confined in a multiply connected region by making use of a spin- 
degree of freedom of the order parameter of a BEC. We demonstrate that this persistent current 
is easily transformed into a vortex. Relaxation processes of these BEC after the confining field is 
turned off are also studied so that our analyses are compared with time of flight experiments. The 
results are shown to clearly reflect the existence of a persistent current. 
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I. INTRODUCTION 

Since the discovery of Bose-Einstein condensation 
(BEC) in alkali atoms 0, numerous attempts are made 
to show that the system exhibits superfluidity. One of 
these attempts is to create and observe quantized vor- 
tices by confining the system in a rotating anisotropic 
trap. Recently the vortex is created ||] using the two- 
component BEC. 

In the present paper, we propose a novel method to cre- 
ate a state with a persistent current or a vortex, where 
the hyperfine spin F of alkali atoms is fully utilized. This 
method has been briefly reported in Q, and more de- 
tailed theoretical and numerical analyese are made in the 
present paper. Although we restrict ourselves mainly to 
the case F = 1 to simplify our discussions, our method 
is also applicable to cases with an arbitrary F. 

BEC with F — I may be expressed in terms of a three- 
component order parameter, similarly to the spin or the 
orbital part of superfluid 3 He. In particular, a spin- 
polarized BEC has the same order parameter as that of 
3 He-A [|. In the case where the spin-exchange interac- 
tion is ferromagnetic, the BEC is spin-polarized even in 
the absence of an external magnetic field. Even when the 
spin-exchange interaction is antiferromagnetic, it may be 
also spin-polarized provided that the Zeeman energy is 
larger than the spin-exchange energy. Accordingly, each 
of the weak field seeking state and the strong field seeking 
state has the same order parameter as in 3 He-A. 

In contrast with 3 He-A, the local order parameter con- 
figuration, known as a texture in 3 He ||, of a spin- 
polarized BEC may be easily controllable by a magnetic 
field. Making use of this property, a BEC in a vortex 
state or with a persistent current can be continuously 
created from a BEC without circulation by adiabati- 



cally changing an external magnetic field as shown below. 
The general theoretical framework for describing a spinor 
BEC Q was given by Ohmi and Machida |(| and indepen- 
dently by Ho [Q, which turned out to be equivalent. This 
framework is based on the Bogoliubov theory extended to 
a vectorial order parameter with three components, cor- 
responding to rrif = 1, 0, — 1 of the F — 1 atomic hyper- 
fine state. As a result, the generalized Gross-Pitaevskii 
(GP) equation has been constructed. They calculated 
low-lying collective modes such as sound wave, spin wave, 
and their coupled modes and predicted various topologi- 
cal defects or spin textures. 

This paper is organized as follows. In the next sec- 
tion, the order parameter of a BEC with F = 1 is dis- 
cussed. We employ two sets of basis vectors and their 
transformations are also considered. Then the general- 
ized Gross-Pitaevskii equation is introduced. In Section 
III, we consider the cross disgyration texture which is ex- 
pected to appear in an Ioffe-Pritchard trap. Section IV is 
the main part of the present paper. We first consider an 
axially symmetric BEC without a current confined in an 
Ioffe-Pritchard trap with an optical plug. A strong mag- 
netic field is applied along the axis of the BEC. Then the 
sign of this axial magnetic field is adiabatically changed 
so that the local magnetization vector flips in the end of 
this process. Then it is shown that a persistent current 
with two units of circulation is created. If the optical plug 
may be turned off at this stage, we are left with a vortex 
line. The topological justification of this behavior is also 
given. Observational consequences of the existence of a 
quantized vortex or a persistent current are discussed in 
Section V, where the relaxation of the order parameter 
after the confining fields are turned off is studied. Section 
VI is devoted to summary and discussions. 
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II. SPINOR BOSE-EINSTEIN CONDENSATE 

A. Spinor order parameter 

Let F a (a — x, y, z) be the angular momentum opera- 
tors with F = 1. The eigenvalues of F z are 1, 0, —1 and 
their corresponding eigenvectors, that satisfy F z \i) = i\i), 
are 



|1> 



|0) 



1) 



(1) 



In this basis, called the ZQ basis, F a are represented as 



F x = 



F, = 



F z = 




(2) 



which satisfy the commutation relation [F a , Fp] = 
iF^e a pj. The order parameter |$) is expanded in terms 
of \i) as 



I*) 



E *< 



(3) 



4 = 0. ±1 



It is convenient for our purposes to introduce another 
set of basis vectors \a) (a = x, y, z) called the XY Z basis 
which satisfy F a \a) = 0. Note that \z) = |0) by defini- 
tion. The vectors \x) and \y) are obtained by rotating \z) 
around the y axis and the x axis by ±7r/2; 



\x) = exp (-i^F y ) \z) = ^(-|1) + I - 1)), 



exp \z) = -L(|l) + |-1)). 



(4) 
(5) 



Then the order parameter |$) may be decomposed in 
terms of la) as 



i*>= E *«i 



(6) 



a—x,y ,z 



The components Vti and are related with each other 
as 




— — 

1 

— — n 

V2 V2 




(7) 



In a spin-polarized BEC, the weak or strong field seek- 
ing state is represented by an order parameter: 



* = 4=e <a (w + m) 
v2 



(8) 



in the XYZ basis, where ip — 

are real unit vectors. We also define 

I = rh x n, 



^i-l 2 and m and h 



(9) 



which represents the direction of the atomic hyperspin. 
The three real vectors (rh,n, l) form a triad. 

In the above explanation the direction of the axis of 
quantization is named z. We may take this direction 
arbitrary. When the axis is parallel to the direction of 
the magnetic field (B in this paper), the Zeeman energy 
term is written most simply. We call this B-quantized 
(BQ) notation. When the direction of the axis does not 
vary spatially and parallel to the z axis, we call this z- 
quantized (ZQ) notation. The kinetic energy term is 
written simply in this way. The numerical details are 
given in Appendix [a]. 



B. Gross-Pitaevskii equations 

The time-dependent form of the Gross-Pitaevskii (GP) 
equation with a spin-degree of freedom obtained by Ohmi 
and Machida |^|, originally in the XYZ basis, can be 
rewritten in the ZQ notation as 

d 



—^ ] ={h jk +g n 5 ]k Y J \^l 



dt 

where 

h jk {r) 



-9s E E (VMlpVp) (F a ) jk }*k (10) 

a Ip 



B = 



h z v 2 

2m 
B z 





d + V(r))6 jk -B jk , 



B x -iB v 



B w +iB y 
V2 





—B. 



(11) 



(12) 



m is the mass of an atom and B = (B x , B y , B z ) is a mag- 
netic field scaled so that the amplitude is the Zeeman en- 
ergy of the atom. The potential V(r) is spin-independent 
and i,j = 0, ±1 are the spin indices in the ZQ basis. The 
parameters g n and g s denote the strength of the inter- 
actions. The relationship between (B x , B y , B z ) and B is 
explained in Appendix Time-independent solutions of 
GP equation are obtained by solving 

= {h jk + g„S 3kz 2\^i\ 2 
i 

+9* E E (*j(^a)lp*p) (F a ) jk }*k. (13) 

a lp 
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The above equations are derived from the Hamiltonian 

3 

+TE*5( r )*t( r )**(r)*i(r) 



jk 



■|e(e«!w(w*w) 

a \ jk / 



(14) 



n 2 v 2 

2m 



M + ^(r) T S(r) + ( ff „+ ffs )|*| 2 |*. 

(19) 



We have ignored the small corrections, which comes from 
the spatial derivative, to the kinetic energy term. This is 
the usual GP equation without the spin degrees of free- 
dom. 

When B z = 0, the order parameter with the highest 
eigenvalue corresponds to the strong field seeking state is 



III. STRONG AND WEAK FIELD SEEKING 
STATES IN IOFFE-PRITCHARD TRAP 

We consider a system of BEC which is uniform along 
the z axis. The cylindrical coordinates r = (r, 4>, z) are 
introduced. Suppose that an Ioffe-Pritchard field 



B = (Bx(r)cos0,-Bj.(r) sm(/>,B z ) 



(15) 



is applied to the system. We treat two-dimensional sys- 
tem (uniform along the z axis) in the following calcu- 
lations and B z is treated as a constant, which differs 
from the usual Ioffe-Pritchard field. There should be a 
blue-detuned laser beam penetrating along the z axis to 
prevent the atoms from escaping from the trap by spin- 
flipping. Ohmi and Machida 0j have shown that there 
appears the cross disgyration when B z = 0. 

Let us derive the configuration of the condensate in 
this system. The 23-matrix Eq. (112) becomes 



( B * B ^ \ 
B 



V2 



V 



B x ^t -B z 



(16) 



T he eigenv alues of B are ±B and where B = 
\J B'± + .Bf , and the corresponding eigenvectors, denoted 
by |^)bq, are 



± 1 




|0>! 



V2B 



(B ± B z )e l4 

±y/2B ± 
(B T B z )e-' 

-Bxe^ 
V2B Z 



(18) 



The vectors |1)bq and | — 1)bq are identified with the 
strong field seeking state and the weak field seeking 
state respectively. Accordingly when the whole system 
is in the strong or the weak field seeking state, the or- 
der parameter is written in terms of these vectors as 
|*) = (r)| ± 1) B q. The £4- term in Eqs. © and @) 
is then simplified to ±B(r)^ so that the GP equation 
takes the form 




1> 



l i<j) 

1 
73 

2 e 



(20) 



where w is an integer and t/j is the amplitude of St. In 
the XYZ basis, this is rewritten as 







(i;) 









-1-7= sm c 
-^cos< 

75 



(21) 



The corresponding rh,n and I vectors are 

rh = (sin(f>sm(w(f>),cos cj)sm(w(f)),co8(w(f))), 
fi= (— sin<j)cos{w<j)), — cos4>cos(w(f)),sm(w<p)), (22) 
I = (cos (j), — sin 4>,0). 



In the weaking field seeking state, the order parameter 
is written as 




l p i<t> 



-*\ ~h 



(23) 



or 



(17) 



















k P —i<t 
2 e 



■ 1 

-t-j= COS( 

v- 1 



Jw<l> 



(24) 



V2 



The corresponding triad is 

rh = (sin 0sin(ii?0), cos <p sm(w<j)) , — cos(u>0)), 
n = {— sin0cos(w0), — cos <jf>cos(«;0), — sin(io^)), (25) 
I = (— cos 0, sin (j>, 0). 



This i vector field in both (|22j) and ( |25| ) clearly represents 
the cross disgyration. 
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IV. CREATION OF VORTEX 

In the present Section, we propose a simple method to 
create a persistent current (and also a vortex state) in a 
torus-shaped BEC from a state with no persistent cur- 
rent. Then it will be shown that this persistent current 
is easily transformed to a vortex. 

In the following we consider two cases separately. In 
Case I, the condensate is confined optically and the spin 
of each atom points the direction of the magnetic field. 
Namely, atoms are strong field seekers. In Case II, the 
magnetic field is also used to confine the condensate. The 
spin of atoms points are antiparallel to the magnetic field 
and the atoms are weak field seekers. The former case has 
an advantage in its theoretical simplicity. On the other 
hand, the latter case does not require apparatus for op- 
tical confinement except for the repulsive plug around 
r — and can be realizable more easily. 

In the following discussions we consider a BEC of 
Na atoms with F = 1. The mass of the atom m — 
3.81 x 10 _26 kg, the interaction parameter g n = 4:7rh 2 a/m, 
and the scattering length a = 2.75 x 10 _9 m are employed. 
We ignore the other interaction parameter g s . This is 
possible since the whole condensate is assumed to be in 
either strong or weak field seeking state as a whole in the 
following. In those states the interaction terms in the CP 
equation are reduced as shown in Eq. ([l9|). 

The particle density is taken to be around 10 19 m -3 
and the detailed density profile is given in each figure. 
The time span of the persistent current creation process 
T = 30ms is chosen since this is between (the Larmor 
frequency) _1 ~ lm sec and the life time of the conden- 
sate ~ 1 sec. 



A. Case I: Optical confinement 

The external magnetic field B(r, </>, z) takes the form 

(B x ,B y ) — B± (cos 0, - sin <j>) , 

B z — B z0 cos(7r(l — t/T)), (26) 
B± = B' ± rsia{Tt{l-t/T)), 

where t is the time. The factors of the field are taken 
to be B' ± = 200ft, [J//im] and B z0 = 2 x 10 3 /i [J] (note 
that the scaled magnetic field represents the Zeeman en- 
ergy in fact). One finds from Eq. ( |26| ) that B z flips from 
— B z o to B z q so that I also flips in the end of the evo- 
lution. The Larmor frequency is ojl ~ 0.6 x 10 3 [Hz] for 
B ~ 2xl0 3 /i[J], Thus lu l xT ~ 18. The spin-independent 
potential is 

V(r)^^lr 2 + Ue X ^-Q, (27) 

with v = 200 [Hz], U = 1 x 10 4 /i[J] and r = 5 [/mi]. The 
first term of Eq.(27) is the confining potential while the 
second term is the potential produced by the optical plug. 



We obtained the order parameter profile by numeri- 
cal integration of the time dependent GP equation ([h]) . 
The initial state is taken to be the ground state with 
no circulation. The magnetic field changes the direction 
slowly from upward to downward according to Eq. ( |2^ ) 
as shown in Fig. |l|(a) so that the atoms remain in the 
strong field seeking state. The change in the number of 
the fc-th component 

N k (t) = [ \V k (v,t)\ 2 d 2 r (A; = -1,0,1) (28) 



is shown in Fig. 0(b). The total particle density 

n(r,t)=Y / \*k(r,t)\ 2 (29) 

k 

changes with the magnetic field as shown in Fig. 0(c). 

The resulting triad configurations are shown in Figs. 
^ for t = 0,15,30 [ms]. Figure ||(a) shows the initial 
vector configurations. Figure ^|(b) shows the vector con- 
figurations when t = 15ms. The Z-texture is nothing but 
the cross disgyration explained in the previous Section 
since B z vanishes now. We see that the vectors m and 
h rotates around I by 27r as we go around the z axis 
once. Finally when t = 30ms, we obtain a texture with Z 
almost points up everywhere. The vectors m and h ro- 
tates around Z by 4tt as one goes around the z axis once 
in this case and one finally obtained a uniform Z-texture 
with a circulation with the winding number 2. 

Now that a persistent current is created, it is easy to 
transform this into a vortex. The BEC has the k = 1 
component only at t = 30ms. Then there are no atoms 
near the axis r = since the centrifugal force prevents 
the atoms to come close to the axis. Thus one may simply 
turn off the optical plug to obtain a vortex. The details 
are analyzed in next section. 
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FIG. 1. The process of the persistent current creation in 
Case I. (a) Magnetic field at r — as a function of time. 

Because B± oc r, the total magnetic field varies slightly at 
r 7^ lO^im. (b) Particle numbers Nk as a function of time. 
The condensate has component only at t = and 
component only at t — T. (c) Total number density distribu- 
tion. The condensate is almost fixed at around r = 9fim by 
the optical potential V(r). However the change of the total 
magnetic field [see caption (a)] causes the change in the radial 
distribution as shown here. 




(a) 




(c) 

FIG. 2. The triad configurations for t = 0, 15, and 30 ms. 
The arrows denote I while m and ft axe on the disk. The line 
on the disk is rh. Note that there is no need to draw n since 
it is uniquely given by I x rh. When t = and 30 ms, I point 
down and up respectively. When t = 15 ms, I lies almost on 
the xy plane. 



B. Case II: Magnetic confinement 

In Case I, the condensate is confined with spin- 
independent optical trap. Here in Case II, we consider 
the situation where the quadrapole magnetic field Eq. 
(|l5| ) always exists and the additional field B z change from 
a large positive value to a large negative value as shown 
in Fig. ||(a) : 

(B x ,B y ) = B±(cos(/), -sin», 

B z =B z0 (l-2t/T), (30) 
B x = B' ± r. 

Here we take B'^ = AOOh [J/jum], B z0 = 2 x 10 4 h and 
T = 30 [ms]. Contrary to Case I, the atoms are in the 
weak field seeking state and the gradient of B r is respon- 
sible for the confinement of the condensate. The optical 
plug produces a spin-independent potential 
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V(r) = U exp 



2rl 



(31) 



where we take U = 1 x 10 4 /i[J] and ro = 5[/im]. 

The evolution of the order parameter field is analyzed 
numerically and it was found that the order parameter 
configurations in this case is essentially the same as in 
Case I. Figure || (b) shows the temporal evolution of the 
components N^tt) and Fig. ||(c) shows the total density 
profile at various time. 



The persistent current created here may also be trans- 
formed into a vortex. The BEC is mostly made of the 
k = 1 component at t — 30ms and there are only a small 
number of atoms near r = 0. Suppose the optical plug 
is turned off. Then atoms will escape but this process 
should be very slow. Thus one expects that the vortex is 
stable for a considerable period of time. 




FIG. 3. The process of the persistent current creation in 
Case II. (a) Magnetic field at r = lO^im as a function of time. 
Initially B z is positive, in contrast with Case I. (b) Particle 
numbers N k as functions of time. The condensate is mostly 
made of \1/ _ i component at t = while mostly of \&i com- 
ponent at t = T. (c) Total density distribution. Contrary 
to Case I, the condensate is confined by the B± part of the 
magnetic field at larger r. Then the radial change is more 
eminent than that of Case I. 
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C. Mathematical analysis of continuous creation of 
circulation 

It may be surprising that we have continuously created 
a persistent current (a vortex) from a system without 
circulation. Mathematically this is justified by invoking 
homotopy theory. Let us denote a rotation around direc- 
tion n by an angle a by a "vector" ae. This rotation is 
expressed as a rotation matrix 



i?(e, a) = (1 — cos a)hihj + cos aSi 



sin as. 



jk h k . (32) 



(b) 



Since a may be restricted within the region < a < n, 
the set of all the rotations is represented by a ball B 3 
with the radius 7r. Note however that the points ire and 
— 7re corresponds to equivalent rotations. Thus all the 
antipodal points on the surface of the ball must be iden- 
tified. This space B 3 /Z2 is called the three-dimensional 
real projective space, denoted by RP 3 . 

Let us take a "standard" triad (rho,noj'o) shown in 
Fig. H Then an arbitrary triad is obtained by applying a 
certain rotation a uniE to the standard triad. Thus the 
local vector configuration is in one-to-one correspondence 
with a point in RP 3 . 

Consider an order parameter configuration shown in 
Fig. |(a). When one circumnavigates the circle, one finds 
that all the triads along the circle are obtained from the 
standard one by applying no rotations, 

namely a — 0. Thus this circle is mapped to the 
origin of RP 3 . Next consider the triads in Fig. ||(b). 
As one goes along the circle, one finds that the stan- 
dard triad is rotated by an angle ir/2 around the axis 



G 



e = (— sin 4>, — cos4>, 0) shown in Fig. 4 Thus this circle 
is mapped to a circle in RP 3 with the radius 7r/2, see 
Fig. 4. Finally, consider the triads in Fig. ||(c). All the 
I- vectors point up and the standard triad is rotated by 
it around the axis e given above. Thus the the circle in 
Fig. 2(c) is mapped to a great circle with the radius it. 

The change of the images in RP 3 , namely a point — * a 
circle with the radius 7r/2 — > a circle with the radius n, is 
continuous (or more precisely homotopic), which shows 
that the deformation of the triads is indeed continuous. 

In the configuration (a), the vector fh does not ro- 
tate around I and therefore there is no current flowing 
around the circle. In configuration (b), however, rh ro- 
tates around I once while one goes along the circle once, 
which implies a vortex of winding number 1. Similarly, 
rh winds around I twice in the configuration (c) and the 
vortex has the winding number 2. We stress again that 
the creation of this winding number (or circulation) is 
continuous and the final configuration is stable so far as 
the external magnetic field forces I to point upward. 




FIG. 4. The real space configuration of e, the RP 3 space 
and the standard triad. The thick arrows are ae vectors. 



D. Angular momentum analysis of current creation 

The continuous current creation can be explained more 
generally using the eigenvectors of the ,6-matrix with the 
highest and lowest eigenvalues in the BQ basis. The mag- 
netic field expressed in S-matrix is the dominant factor 
to determine the behavior of the order parameter '4'. Wc 
first discuss the BEC with F = 1, which we have ana- 
lyzed in the present paper. 

The magnetic fields Eq. ( ^6|) i n Case I and Eq. ([30|) in 
Case II are of the form Eq. (|15|). The corresponding B- 



matrix is given by Eq. (Jlq). Because the magnetic fields 
are sufficiently strong, we may assume that the order pa- 
rameter \I/ is proportional to the highest (upper sign) or 
lowest (lower sign) eigenvector 

/ (B±B z )e^ \ 
| ± 1) BQ = ^ ±V2B ± . (33) 
ZB V (B T B z )e-* ) 

of the K-matrix. Case I uses the highest (upper sign) 
and Case II uses the lowest (lower sign) eigenvector. The 
order parameter is 

\ / (B±B z y* \ 

*o = C ±V2B ± e lw t (34) 

*-i / \(BTB z )e-^J 

where C is a complex number independent of cj> and w is 
an integer. 

Let us consider Case I. (Case II can be analyzed sim- 
ilarly.) When t = 0, the magnetic field is (B±,B Z ) = 
(0, tB z q). As shown in Eq. (£34|), the condensate behaves 
as *_i cx e^- 1 ^ and * = *i = 0. Because we start 
with a state with no circulation, the integer w must be 
1. 

The components B z and B± change as shown in, for 
example, Fig. [j](a) from t = to t = T. Both B± and 
(B ± B z ) have finite values during the change and the 
condensate stays in the strong field seeking state. The 
phase factor w will not change during the process for the 
order parameter to be defined uniquely. Thus we take 
w = 1 throughout the process. When t — T, the mag- 
netic field is (B±,B X ) = (0,±B z0 ). Since w = 1, Eq. 
( |34| ) leads to the conclusion that we have a condensate 
with *i oc e 2 ^ and * = = 0. 

Accordingly a vortex with the winding number 2 has 
been created. 

A similar discussion is possible in the system wit h 
F = 2 atoms using the eigenvector given in Eq. (|A7j). 
Starting from a state with no winding number, we even- 
tually obtain a state with the winding number 4. 



V. DETECTION OF VORTEX: TIME OF FLIGHT 
IMAGING 

The detection of a vortex (or persistent current) has 
been a problem as difficult as their creation. We consider 
the relaxation of the spinor texture after the confining 
field and the optical plug are turned off to facilitate the 
comparison between our theory and experiments, in par- 
ticular the time of flight (TOF) analysis. 

The temporal evolution of the BEC is described by the 
time-dependent GP equation (10). We consider Case I 
and Case II separately. 
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A. Case I 



We consider three cases where the confining potentials 
are turned off separately at t — 0, T/2 and T ms. There 
is no vortex at t = while there is a vortex with the 
winding number 2 at t = T. The comparison between 
these two cases is essential to observe our vortex. 

(i) t — 0: The condensate has a component only. 
The density profile at t = is determined by solving the 
time-independent GP equation. The relaxation process 
after the potentials are turned off is found by solving the 
time-dependent GP equation Eq. (|Io|), whose result is 
shown in Fig. ||(a). Since has no singularity at the 
origin, the condensate fills the central region (r ~ 0) in 
later time, ft is interesting to note that the components 

and $i do not appear in later time since the total 
spin must be conserved. 

(ii) t = T/2: The cross disgyration appears in this 
stage. The order parameter of this texture is given by 
Eq. ( pfj| ) with w = 1. Thus all the components are non- 
vanishing in this case. After the potentials are turned off 
at t = 0, the order parameter relaxes as shown in Fig. 
^|(b). The component has winding number 1 while 
'i'l has winding number 2, and hence they cannot fill the 
central region. The central region near r — may be 
filled only with 'J'-i component since it has vanishing 
winding number. Note also that the component is 
dominant in the vicinity of r ~ f/im. 

(iii) t = T: The vector I points up now and hence 
= * = while \&i 7^ 0. Figure |(c) shows the 

temporal evolution of the order parameter after the po- 
tentials are turned off. Since the order parameter has 
a nontrivial phase factor, it cannot fill the central re- 
gion at all in later times. Similarly to the case (i), the 
components and ^-i do not appear in the relaxation 
process. The absence of the condensate at r — at an 
arbitrary time is a clear distinction between the case (iii) 
and the rest, which may be used to show the existence of 
the vortex or the persistent current experimentally. 

The vacuum region around r — which shows the 
existence of the vortex have length scale around 10/im 
after 7ms relaxation, and the length will be sufficient 
to observe experimentally. Because the length scale of 
the central vacuum region is almost same as that of the 
density waves at larger r (outer), the resolution of the 
imaging will be checked by the outer density waves. 



FIG. 5. Temporal evolution of the condensate, (a) The po- 
tentials are turned off at t = 0. The solid line is the total 
density purely made of |*-i| 2 . (b) The potentials are turned 
off at t = T/2. The solid line is the total density J2 k l* fc | 2 
while the dashed line, the dotted line, and the dashed dot- 
ted line are components l^il 2 , l^ol an d |^-i] 2 ; respectively. 
The corresponding winding numbers are 2, 1 and 0. (c) The 
potentials are turned off at t = T. The solid line is the total 
number density made of |^i| 2 only. The condensate has the 
winding number 2 and cannot fill the central region. This 
should be compared with (a) and (b). 

B. Case II 

Figure ^(a) shows the relaxation process when the po- 
tentials are turned off at t = while they are turned 
off at t = T in Fig. ^|(b). Because the magnetic field is 
not exactly parallel or antiparallel to the z axis, the non- 
dominant component of the condensate appears slightly 
in both cases. 




(a) 



(b) 

FIG. 6. Temporal evolution of the condensates when the 
potentials are turned off at (a) t = and (b) t — T. In both 
cases the description of the lines are the same as that in Fig. 5. 

VI. SUMMARY AND DISCUSSIONS 

In summary, we have proposed a new method to cre- 
ate a persistent current and a vortex with the winding 
number 2 in a Bose-Einstein condensate of alkali atoms. 
The dynamics of vortex creation are simulated by solv- 
ing the time-dependent Gross-Pitaevskii (GP) equation. 
The continuity of this process is justified by invoking ho- 
motopy theory and also by the angular momentum anal- 
ysis. The existence of the vortex may be demonstrated 
by comparing the time of flight (TOF) data before and 
after the vortex creation. 

It is also possible to create a persistent current or a 
vortex with the winding number 1. Suppose one pre- 
pares the ground state in the Ioffe-Pritchard field. The 
resulting texture is the cross disgyration with no wind- 
ing of the rh- vector around the 1- vector as shown in ||. 
Then apply a strong magentic field B z either parallel to 
or antiparallel to the z axis. The Z-vector in the resulting 
texture points up or down, depending on the direction of 
B z or whether the state is weak field seeking or strong 
field seeking. In any case, the m-vector rotates around 
I once as one circumnavigates around the z axis once. 
Thus one obtains a persistent current or a vortex of the 
winding number 1 by simply preparing a sample in the 
Ioffe-Pritchard trap and applying a strong magnetic field 
along the z axis. 

The work of MN is supported partially by the Grant- 
in- Aid for Scientific Research Fund of the Ministry of 
Education, Science, Sports and Culture, No. 11640361. 




APPENDIX A: MAGNETIC FIELD MATRIX 



We consider a system of atoms with spin F = f. 

The F-matrices, which are the angular momentum operators, are given by 
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(F z )jj = f + l-j (j = l,2,...,2/ + l). 

They satisfy the commutation relation [F a ,F^] = iF 7 e a| g 7 . 
We write the magnetic field with the -B vector 
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where < 9 y < ir and < 9 Z < 2ir. The amplitude \B\ is scaled so that it represents the Zeeman energy. 

The order parameter of a BEC is written with a vector of 2/ + 1 components and operators are expressed in terms 
of (2/ + 1) x (2/ + 1) square matrices. We call the operator for the Zeeman energy B- matrix. We can choose the 
quantization axis so that the 2?-matrix is proportional to the matrix F z . We call this the B quantized (BQ) notation 
because B is proportional to the quantization axis. 

The B matrix in the ZQ notation is obtained by successive spatial rotations 9 y along the y axis and 9 Z along the z 
axis as 



B ZQ = U*BbqU, 
U* = exp{-iF z 6 z )ex.p(-iF y v ). 



Let us study a few examples. When / = 1, 
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The eigenvector with the lowest eigenvalue is 
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When / = 2, 
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The eigenvector with the lowest eigenvalue is 
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